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We have observed the rare decay B° — > p + p~ in a sample of 89 million BB pairs recorded with 
the BABAR detector. The number of observed events is 88^21 ±9, with a significance of 5.1 standard 
deviations with systematic uncertainties included. The branching fraction and the longitudinal 
polarization are measured to be B(B° -* p + p~) = (25±£tf) x 10" 6 and T L /T = 0.98to;°g ± 0.03, 
respectively. 

PACS numbers: 13.25.Hw, 11.30.Er, 12.15.Hh 



Charmless P-meson decays provide an opportunity to 
measure the angles of the unitary triangles constructed 
from the elements of the Cabibbo-Kobayashi-Maskawa 
(CKM) quark- mixing matrix [lj . There has been interest 
in the study of B — > 7T7t and pit decays, where the time- 
dependent CP-violating asymmetries arc related to the 
CKM angle a = arg [ ~V td V t * b / V ud V* b ] , and interference 
between tree and loop (penguin) amplitudes could give 
rise to direct CP violation. The decay B Q — ► p + p~ is 
another promising mode for CP- violation studies and has 
the advantage of a larger expected decay rate and smaller 
uncertainty in penguin contributions. The measurements 
of the amplitudes in B decays to two vector particles 
provide additional tests of theoretical calculations 0, 0, 

The decay B° — > p + p~ is expected to proceed through 
the tree-level b — > u transition and through CKM- 
suppressed b — > d penguin transitions, as illustrated in 
Fig. H 0- The extraction of a from measurements 
made with this decay requires an understanding of the 
contributing amplitudes. It also requires proper account- 
ing for CP-even (S- and D-wave) and CP-odd (P-wave) 
components in the decay amplitude. The recent limit on 
the B° — > p°p° decay rate |6j] and the measurements of 
the B + — ► p + p° branching fraction 0, Q place experi- 
mental limits on the contribution of penguin amplitudes. 
Measurements of the longitudinal polarization, defined as 
the ratio between the longitudinal and total decay rates 
/l = r^/r 2], in the B + — > p + p a decay provide evidence 
that the CP-even component dominates in B —* pp de- 
cays 0,0. 

In this paper we report the observation of the B a — > 
p + p~ decay mode and measurements of its branching 
fraction and the amount of longitudinal polarization in 
the decay. We also make a quantitative estimate of pen- 
guin contributions in this decay using our earlier mea- 
surements in isospin-related B — ► pp modes. 

We use data collected with the BABAR detector Q at 
the PEP-II asymmetric-energy e + e~ storage ring. These 
data represent an integrated luminosity of 81.9 fb _1 at 
the e + e~ center-of-mass (CM) energy of the T(4S') res- 
onance (y/s = 10.58 GeV, on-resonance) , corresponding 
to 88.9 million BB pairs, and 9.6 fb _1 approximately 
40 MeV below this energy (off-resonance). 

Charged-particle momenta are measured in a track- 
ing system consisting of a five-layer double-sided silicon 
vertex tracker (SVT) and a 40-layer central drift cham- 
ber (DCH), both situated in a 1.5-T axial magnetic field. 




d d d d 



(a) (b) 

FIG. 1: Diagrams describing the decay B° — + p + p~ ': (a) 
dominant tree diagram, (b) gluonic penguin diagram. 



BABAR achieves an impact parameter resolution of about 
40 jirsx for the high- momentum charged particles from the 
B decay, allowing the precise determination of decay ver- 
tices. The tracking system covers 92% of the solid angle 
in the CM frame. 

Charged-particle identification is provided by measure- 
ments of energy loss (dE/dx) in the tracking devices 
(SVT and DCH) and by an internally reflecting ring- 
imaging Cherenkov detector (DIRC). A K-tt separation 
of better than four standard deviations (a) is achieved for 
momenta below 3 GeV, decreasing to 2.5 a at the highest 
momenta in the B decay final states. Photons are de- 
tected by a CsI(Tl) electromagnetic calorimeter (EMC). 
The EMC provides good energy and angular resolution 
for detection of photons with energy in the range 20 MeV 
to 4 GeV. The energy and angular resolutions are 3% and 
4 mrad, respectively, for a 1 GeV photon. 

Hadronic events are selected based on track multiplic- 
ity and event topology. We fully reconstruct B° — > p + p~ 
candidates from the decay products of the p ± — > 7r ± 7r° 
and 7T° — > 77 decays. Charged-track candidates are 
required to originate from the interaction point, have 
at least 12 DCH hits and have a minimum transverse 
momentum of 0.1 GeV. Charged-pion tracks are distin- 
guished from kaon and proton tracks with a likelihood 
ratio that includes dE/dx information from the SVT and 
DCH, and, for momenta above 0.7 GeV, the Cherenkov 
angle and number of photons measured by the DIRC. 
Charged pions are distinguished from electrons primar- 
ily on the basis of their EMC shower energy and spatial 
profile. 

We reconstruct tt° mesons from pairs of photons. Pho- 
ton candidates are required to have a minimum energy of 
30 MeV, have a shower shape consistent with the photon 
hypothesis, and not be matched to a track. The typi- 
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FIG. 2: Definition of helicity angles Ox, 82, and <f>, for the 
decay B° —* p + p~ . The tt ± ty° final states are shown in the 
p rest frames. 



cal experimental resolution for the measured tt° mass is 
7 MeV. We require ir° candidates to have an invariant 
mass within 15 MeV of the true tt° mass. The invariant 
mass of the p^ candidate (m^± n o) is required to be in 
the range 0.52 to 1.02 GeV. The helicity angles 8\ and 
62 of p + and p~ are defined as the angles between the 
7r° direction and the direction opposite the B in each p 
rest frame as shown in Fig. [5] The helicity angles are 
restricted to the region —0.75 < cos 812 < 0.95 to sup- 
press combinatorial background and reduce acceptance 
uncertainties due to low-momentum pion reconstruction. 

The B meson candidates are identified from two nearly 
independent kinematic observables ||, the beam energy- 
substituted mass m E s = [(s/2 + Pi ■ PB) 2 /Ef — pj] 1 / 2 
and the energy difference AE = (EiEB — Pi ■ Pb — 
s/2)/y/s, where (Ei,Pi) is the e + e~ initial state four- 
momentum, and (Eb,Pb) is the four-momentum of the 
reconstructed B candidate, all defined in the labora- 
tory frame. For signal events, the mgs distribution 
peaks at the B mass and the AE distribution peaks 
near zero. Our selection requires toes > 5.2 GeV and 
\AE\ < 0.2 GeV, while the signal resolution is roughly 
3 MeV and 50 MeV, respectively. The sideband re- 
gions are defined as 5.2 GeV < toes < 5.27 GeV or 
0.1 GeV < \AE\ < 0.2 GeV. 

To reject the dominant continuum background (from 
e + e~ — > qq events, q = u,d,s,c), we require |cos#t| < 
0.8, where 8t is the angle between the thrust axis of the 
B candidate and the thrust axis of the rest of the tracks 
and photon candidates in the event, calculated in the CM 
frame. The distribution of | cos 0t\ is sharply peaked near 
1.0 for jet-like events originating from qq pairs and nearly 
uniform for the isotropic decays of the B meson. A Fisher 
discriminant (T) combines 11 observables: the polar an- 
gle of the B momentum vector and the polar angle of the 
.B-candidate thrust axis, both calculated with respect to 
the beam axis in the CM frame, and the scalar sum of the 
CM momenta of charged particles and photons (exclud- 
ing particles from the B candidate) entering nine coaxial 
angular intervals of 10° around the i3-candidate thrust 
axis 0- 

The selected sample contains 54,042 events most of 



which populate sidebands of the observables. Back- 
ground from other B decays is estimated with Monte 
Carlo (MC) simulation [l£(; it contributes 5% of the 
events in the selected sample. This background compo- 
nent, arising mainly from b — > c transitions, is explicitly 
included in the fit described below. 

We use an unbinned, extended maximum-likelihood fit 
to extract simultaneously the signal yield and polariza- 
tion. There are three event categories j: signal, contin- 
uum qq, and BB combinatorial background. The likeli- 
hood for each B° — > p + p~ candidate i is defined as 

3 

C i = Y,n j V j {SvJ), (1) 

3=1 

where each of the Vj{xi](3) is the probability density 
function (PDF) for seven observables afj (toes, AE, T , 
m^+^o, m n - 7r o, 0i, 62) and is described by the PDF pa- 
rameters p. The event yields tij for each category j are 
free parameters in the fit. We allow for multiple can- 
didates in a given event by assigning to each selected 
candidate a weight of 1/iVj , where Ni is the number of 
candidates in that event. The average number of candi- 
dates per event is 1.27. MC simulation shows that this 
procedure does not introduce bias while providing a small 
statistical improvement over the random choice of a can- 
didate in a given event. The extended likelihood for a 
sample of iV can d candidates is 

C = exp L- J2 n jj f[ ex P (^r ) • ( 2 ) 

The correlations among the input observables 
found to be small for both the background (<5%) and 
signal (<10%), except for angular correlations in the sig- 
nal. The Vj(xi, (3), for a given candidate i, is the product 
of PDFs for each of the observables and a joint PDF for 
the helicity angles, which accounts for the angular corre- 
lations in the signal and for detector acceptance effects. 
We integrate over the angle </> between the two decay 
planes shown in Fig.|2 leaving a PDF that depends only 
on 6\, 62, and the unknown longitudinal polarization Jl- 
The differential decay rate 0] is 

I d 2 r _ 

r d cos 6*i d cos 82 
/l) sm 2 #1 sin 2 2 + / L cos 2 ^ cos 2 2 1 . (3) 

The PDF parameters (3, except for f^, are ex- 
tracted from MC simulation and on-rcsonance Toes & n d 
AE sidebands, and are fixed in the fit. The res- 
olutions are adjusted by comparing data and simu- 
lation in calibration channels with similar kinematics 
and topology, such as B — > Dp + ,Dn + with D — > 
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K+ir- (ir°), K°ir- {ir°), K+ir-ir- , K°ir-ir+ . To describe 
the signal distributions, we use Gaussian functions for 
the parameterization of the PDFs for rngs and AE, and 
a relativistic P-wave Breit-Wigner distribution for the 
p ± resonance masses. The angular acceptance effects are 
parameterized with empirical polynomial functions for 
each helicity angle and are included in the joint helicity- 
anglc PDF as a factor multiplying the ideal distribution 
in Eq. ©. 

For the background PDFs, we use polynomials or, in 
the case of tties, an empirical phase-space function 
In the background PDF we incorporate a small linear 
correlation between the curvature £ of the phase-space 
function and the value of T . The background param- 
eterizations for the p ± candidate masses also include a 
resonant component to account for p production. The 
background helicity-angle distribution is also separated 
into contributions from combinatorial background and 
from real p mesons, both described by polynomials. For 
both signal and background, the PDF for T is repre- 
sented by a Gaussian distribution with different widths 
above and below the peak. 

PDF parameters for the background from other B de- 
cays are determined from MC simulation. The contri- 
bution from charmless B decays with similar topology 
(cross-feed modes) such as B — * pir, p°p + , pK* , ain, 
and aip is estimated with MC modeling and is fixed in 
the fit. Each branching fraction for the cross-feed modes 
is estimated to be in the range (1-3) x 10~ 5 . The branch- 
ing fractions for these and many other modes are taken 
from the most recent measurements 0, 0, or extrap- 
olated from other results with a flavor-SU(3)-symmetry 
approximation. 

The selected B° — > p + p~ events fall into three cate- 
gories. MC simulation of events with longitudinal po- 
larization shows that roughly 30% of the events contain 
only misreconstructed candidates. Approximately 20% 
of the events contain both correctly and incorrectly re- 
constructed candidates. The remainder contain only cor- 
rect candidates. Misreconstruction occurs when at least 
one candidate photon in a tt° candidate or one charged 
track in a p candidate belongs to the decay products of 
the other B. The distributions that show peaks for cor- 
rectly reconstructed events have substantial tails, with 
large uncertainties in MC simulation, when misrecon- 
structed events are included. These tails would reduce 
the power of the distributions to discriminate between 
the background and the collection of correctly and incor- 
rectly reconstructed events. We choose, therefore, to rep- 
resent only the correctly reconstructed candidates in the 
signal PDF. Misreconstructed candidates are predomi- 
nantly accommodated by the combinatorial background 
PDF. Fitting to determine the number of correctly re- 
constructed candidates has an efficiency less than 100% 
since some fraction of the events have both correctly and 
incorrectly reconstructed candidates. Monte Carlo simu- 



TABLE I: Summary of the fit results; n s i g is the fitted number 
of signal events, S is the significance, /x, is the longitudinal 
polarization, e denotes the reconstruction efficiency, and B 
is the branching fraction of the B° — > p + p~ decays. The 
first uncertainty is statistical and the second systematic. The 
efficiency (e) and significance (5) include systematic uncer- 
tainties, and the significance without systematics is given in 
parentheses. 



Quantity Measured Value 

n sig 88±2? ± 9 

S 5.1 cr (5.5 <t) 

H 0.98 tom ± 0.03 

3.9t8;S % 

B (25 ±1 ±1) x 10-° 



lation finds this efficiency to be 87%. 

In this analysis, we do not include a fit component 
for other B decays with the same final-state particles se- 
lected within the p resonance mass window, such as non- 
resonant decays B° — > 7r + 7r~7r°7r° and B° — > p ± 7r T 7r°. 
The contribution of these decays to the fit results is sig- 
nificantly suppressed by the selection requirements on the 
masses and by the mass and helicity-angle information 
in the fit; they are examined in the context of mass and 
helicity-angle distributions, as discussed below. 

The event yields rij and polarization fa are obtained 
by minimizing the quantity \ 2 = — 21n£. The depen- 
dence of x 2 on a fit parameter rij or fa is obtained with 
the other fit parameters floating. Their values are con- 
strained to the physical range rij > and < fa < 1. 
Statistical uncertainties correspond to a unit increase in 
X 2 ■ The statistical significance of the signal is defined as 
the square root of the change in x 2 when the number of 
signal events is constrained to zero in the likelihood fit. 

The results of our maximum-likelihood fits are sum- 
marized in Table [i] The statistical significance of the 
B° — > p + p~ signal is 5.5 a. We find that the p ± mesons 
in B° — > p + p~ decays are almost fully longitudinally po- 
larized. To compute the branching fraction, equal pro- 
duction rates for B°B° and B + B~ are assumed. To check 
the stability of our results we refit, removing each observ- 
able from the fit in turn, and find consistent results. The 
measured uncertainties in the number of fitted events and 
the polarization, the statistical significance, and the fit 
X 2 value are well reproduced with generated MC samples. 

The projections of the fit input observables are shown 
in Fig. |2| The projections are made after a re- 
quirement on the signal-to-background probability ratio 
V s i s (xi; P)/Vbk B {xi; 0), where V s i g and 7\k g are the sig- 
nal and the dominant continuum background PDFs de- 
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FIG. 3: Projections onto the observables ities, AE, m^+^o, 
m^-^o, cos Si, and cos 62 after a requirement on the signal- 
to-background probability ratio V a i s /Vbk s with the PDF for 
the plotted observable excluded. The points with error bars 
show the data, the solid (dashed) line shows the signal-plus- 
background (background only) PDF projection. 



fined in Eq. , but with the PDF for the plotted observ- 
able excluded. The points with error bars show the data 
with (40-60)% of signal retained, while the lines show 
the corresponding PDF projections. 

To check the sensitivity of our results to the presence 
of nonrcsonant B° — > 7r + 7r~7r°7r° and B° — ► p ± 7r =F 7r° 
decays, we explicitly include a fit component for them, 
assuming a phase-space decay model. The selection re- 
quirements alone suppress the B — > 47r (B — > pirn) ef- 
ficiency by two (one) orders of magnitude relative to 
B° — > p + p~ ■ The fit results with a nonresonant com- 
ponent indicate a potential B — > pixix contribution of 
(10 ± 10)% (statistical uncertainty only) of our nominal 
B° — ► p + p~ event yield in Table HJ interference effects 
between the resonant and nonresonant components were 
ignored in this fit. The hypothesis that all the signal is 
nonresonant B — > 4tt (B — > pirir) is excluded with 5.1 a 
(4.4 er) statistical significance. These results are consis- 
tent with our assumption that the nonresonant contribu- 
tion is negligible. 

The systematic uncertainty in the fitted number of sig- 
nal events (n s i g ) originates from the uncertainty in the 
cross-feed £?-decay modeling, which was studied with MC 
generated samples and estimated to be half of the varia- 
tion with cross- feed set to zero (3% uncertainty in n s i g ). 
Systematic uncertainties in the fit originate from assump- 
tions about the PDF parameters. Uncertainties in the 



PDF parameters arise from the limited number of events 
in the background sideband data and signal control sam- 
ples. We vary them within their respective uncertainties, 
and derive the associated systematic uncertainty on the 
event yield (9%). The signal remains statistically signifi- 
cant with these variations (5.1 a including systematics). 

The systematic uncertainties in the efficiency (e) are 
due to track finding (2% for two tracks), particle identi- 
fication (2% for two tracks), and 7r° reconstruction (13% 
for two 7r°s). The fit efficiency is less than 100% because 
of misreconstructed signal events. This has an additional 
systematic uncertainty due to uncertainties in the mod- 
eling of misreconstructed events. We account for this 
with a systematic uncertainty on the efficiency of 7%, 
which is half of the inefficiency; the fit efficiency cannot 
exceed 100% and the frequency of multiple candidate se- 
lection is estimated in the B decay control samples. The 
reconstruction efficiency depends on the decay polariza- 
tion. We calculate the efficiencies using the measured 
polarization and assign a systematic uncertainty ( _J%) 
corresponding to the total polarization measurement un- 
certainty. Smaller systematic uncertainties arise from 
event-selection criteria, MC statistics, and the number 
of produced B mesons. 

For the polarization measurement (/l), we include sys- 
tematic uncertainties from PDF variations that account 
for uncertainties in the detector acceptance, estimated 
with MC, and background parameterizations. This re- 
sults in a total absolute uncertainty of 0.025. The bi- 
ases from the resolution in helicity-angle measurement 
and dilution due to the presence of the misreconstructed 
combinations are studied with MC simulation and give a 
systematic uncertainty of 0.02. 

Observation of the B° — > p + p~ decay completes a 
first set of measurements of the isospin-related B — > pp 
modes 0, Q . The measured branching fraction is consis- 
tent with recent predicted values in the range (18-35) x 
10 -6 Q and the dominant longitudinal polarization im- 
plies a suppression of the transverse amplitude, which is 
expected to be suppressed by a factor of m p /ms p|. The 
rates of the B° — ► p + p~ and B + — ► p°p + decays appear 
to be larger than the corresponding rates of B — > 7T7T de- 
cays 0| • At the same time, the recent measurement of 
the B + — > p°K* + branching fraction does not show 
significant enhancement with respect to B — > irK de- 
cays ^2 , both of which are expected to be dominated by 
b — ► s penguin diagrams. We can use flavor SU(3) to re- 
late b —>■ s and b — ► d penguins analogous to Fig.j^b) [l3| : 
the measured branching fractions indicate that the rel- 
ative penguin contributions in the B — > pp decays are 
smaller than in the B — 7T7t case. 

We make a more quantitative estimate of penguin con- 
tributions in B pp decays using our previous measure- 
ments of B° — » p°p® and B + — > p + p° branching fractions 
and polarization [6j. Since the tree contribution to the 
B° — > p°p° decay is color-suppressed, the decay rate is 
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sensitive to the penguin diagram analogous to Fig. ^b). 
Using the earlier BABAR measurements 0, we obtain a 
90% confidence level (C.L.) upper limit on the ratio of 
the longitudinal amplitudes in the B — > pp decays: 



\A L (B a ^p a p a )\ 2 + \A L (B° 



P°P°) 



2 x \A L (B^ 



B{B°^p a p°) x /l(B° — > p°p°) 



B(B^ 



x f L {B+ -» fflp+) 



< 0.10. 



(4) 



In the above calculation we conservatively assume that 
the B° — > p°p° decay polarization is fully longitudinal 
(/l = 1), use the average branching fraction measure- 
ments for the B and B decays, and assume \Al(B + — > 
= I^l(B~ — > with only a tree-diagram 

contribution. The limit in Eq. Q corresponds to a 
19° uncertainty (at 90% C.L.) on a due to penguin 
contributions in the time-dependent measurements with 



longitudinally-polarized B 



p + p decays, assuming 



isospin relations analogous to those discussed in the con- 
text of B — > 7T7T |14| and ne glec ting the nonresonant and 
1=1 isospin contributions |15| . 

In summary, we have observed the decay B° — ► p + p~ , 
measured its branching fraction B = (25lglg) x 10~ 6 , 
and determined the longitudinal polarization fraction 
/l = 0.98lQ;Qg±0.03. Our quantitative estimates of pen- 
guin contributions in B° — > p + p~ decays and the domi- 
nance of the CP-even longitudinal polarization make this 
decay a promising channel for the measurement of the 
CKM angle a. 

We are grateful for the excellent luminosity and ma- 
chine conditions provided by our PEP-II colleagues, and 
for the substantial dedicated effort from the computing 
organizations that support BABAR. The collaborating 
institutions wish to thank SLAC for its support and 
kind hospitality. This work is supported by DOE and 
NSF (USA), NSERC (Canada), IHEP (China), CEA and 
CNRS-IN2P3 (France), BMBF and DFG (Germany), 
INFN (Italy), FOM (The Netherlands), NFR (Norway), 
MIST (Russia), and PPARC (United Kingdom). Indi- 
viduals have received support from the A. P. Sloan Foun- 



dation, Research Corporation, and Alexander von Hum- 
boldt Foundation. 



[6 

[<: 

[8 
[9 
[10 

[11 

[12 
[13. 
[14 

[is; 



Also with Universita della Basilicata, Potenza, Italy 
Also with IFIC, Instituto de Ffsica Corpuscular, CSIC- 
Universidad de Valencia, Valencia, Spain 
Deceased 

N. Cabibbo, Phys. Rev. Lett. 10, 531 (1963); 
M. Kobayashi and T. Maskawa, Prog. Theor. Phys. 49, 
652 (1973). 

G. Kramer and W.F. Palmer, Phys. Rev. D 45, 193 
(1992). 

R. Aleksan et al, Phys. Lett. B 356, 95 (1995). 
D. Ebert, R.N. Faustov, and V.O. Galkin, Phys. Rev. D 
56, 312 (1997); A. Ali, G. Kramer, and C.-D. Lu, Phys. 
Rev. D 58, 094009 (1998); Y.-H. Chen et al, Phys. Rev. 
D 60, 094014 (1999); H.-Y. Cheng and K.-C. Yang, Phys. 
Lett. B 511, 40 (2001). 

We use the name of one member of a charge conjugate 
pair to refer to both unless explicitly stated otherwise. 
BABAR Collaboration, B. Aubert et al, Phys. Rev. Lett. 

91, 171802 (2003). 

Belle Collaboration, J. Zhang et al, |hep-e x/0306007 
submitted to Phys. Rev. Lett. (2003). 
BABAR Collaboration, B. Aubert et al, Nucl. Instrum. 
Methods A479, 1 (2002). 

CLEO Collaboration, D.M. Asner et al, Phys. Rev. D 
53, 1039 (1996). 

The BABAR detector Monte Carlo simulation is based on 
GEANT4: S. Agostinelli et al, Nucl. Instrum. Methods 
A506, 250 (2003). 

With x = 2mEs/\/s and £ a parameter to be fit, f(x) oc 
xy/1 — x 2 exp [— £(1 — x 2 )] ; see ARGUS Collaboration, 

H. Albrecht et al, Phys. Lett. B 241, 278 (1990). 
Particle Data Group, K. Hagiwara et al., Phys. Rev. D 
66, 010001 (2002). 

M. Gronau and J.L. Rosner, Phys. Rev. D 65, 013004 
(2002). 

M. Gronau and D. London, Phys. Rev. Lett. 65, 3381 
(1990); Y. Grossman and H. Quinn, Phys. Rev. D 58, 
017504 (1998). ^ 

A.F. Falk et al., |hep-ph/0310242| submitted to Phys. 
Rev. D (2003). 



